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Abstract 



fNj ■ The Faddeev equation for three-body scattering below the three-body 

breakup threshold is directly solved without employing a partial wave de- 
£S) \ composition. In the simplest form it is a three-dimensional integral equation 

in four variables. From its solution the scattering amplitude is obtained as 
function of vector Jacobi momenta. Based on Malfliet-Tjon type potentials 
differential and total cross sections are calculated. The numerical stability of 
the algorithm is demonstrated and the properties of the scattering amplitude 
discussed. 



I. INTRODUCTION 



Traditionally three-nucleon scattering calculations are carried out by solving Faddeev 
equations in a partial wave truncated basis. A partial wave decomposition replaces the 
continuous angle variables by discrete orbital angular momentum quantum numbers, and 
thus reduces the number of continuous variables, which have to be discretized in a numerical 
treatment. For low projectile energies the procedure of considering orbital angular momen- 
tum components appears physically justified due to arguments related to the centrifugal 
barrier. However, the algebraic and algorithmic steps to be carried out in a partial wave 
decomposition can be quite involved when solving the Faddeev equations. If one considers 
three-nucleon scattering at a few hundred MeV projectile energy, the number of partial waves 
needed to achieve convergence proliferates, and limitations with respect to computational 
feasibility and accuracy are reached. The amplitudes acquire stronger angular dependence, 
which is already visible in the two-nucleon amplitudes, and their formation by an increas- 
ing number of partial waves not only becomes more tedious but also less informative. The 
method of partial wave decomposition looses its physical transparency, and the direct use 
of angular variables becomes more appealing. 

It appears therefore natural to avoid a partial wave representation completely and work 
directly with vector variables. This is common practice in bound state calculations of few- 
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nucleoli systems based on variational |J and Green's function Monte Carlo (GFMC) methods 
|2|-[5|], which are carried out in configuration space. 

Our aim is to work directly with vector variables in the Faddeev scheme in momentum 
space. In we demonstrated for Malfliet-Tjon type potentials that the two-body Lippmann- 
Schwinger equation can readily be solved in momentum space as function of vector momenta. 
At intermediate energies the strong forward peaking of the i-matrix was quite easily achieved 
through the angular variable but required relatively many partial wave contributions. The 
choice of momentum vectors as adequate variables is also suggested by the nucleon-nucleon 
(NN) force. Here the dependence on momentum vectors can be rather simple, e.g. in the 
widely used one-boson-exchange force, whereas the partial wave representation of this force 
leads to rather complicated expressions 0. 

In || we showed that the bound state Faddeev equation has a rather transparent struc- 
ture when formulated with vector variables compared to the coupled set of two-dimensional 
integral equations obtained in a partial wave decomposed form. Based on Malfliet-Tjon type 
interactions it was demonstrated that the numerical solution of the bound state equation 
using vector variables is straightforward and numerically very accurate. 

In this article we want to show that the solution of the three-body scattering equation 
can also be obtained in a straightforward manner when employing vector variables, i.e. mag- 
nitudes of momenta and angles between the momentum vectors. In this work we concentrate 
on scattering below the three-body breakup threshold. Though we avoid the singularity of 
the free three-body propagator, we already encounter the two-body fragmentation cut re- 
lated to the pole in the two-body f-matrix. As a further simplification we neglect spin and 
iso-spin degrees of freedom and treat three-boson scattering. The interactions employed are 
of Yukawa type, and no separable approximations are involved. The Faddeev equation for 
three identical bosons is solved exactly as function of momentum vectors. To the best of 
our knowledge this is the first time such an approach is carried out. 

This article is organized as follows. Section II describes our choice of momentum and 
angle variables for the unknown amplitude in the Faddeev equation and the integral kernel 
of that equation. The calculation of the only scattering observable, the differential cross 
section is also derived. In Section Ilia we discuss details of our algorithms and numerical 
procedures and present our results. In addition properties of the Faddeev amplitude are 
displayed. In Section IV we discuss an alternative choice of variables and demonstrate that 
this choice leads to the same result as the choice of Section II. We conclude in Section V. 

II. THREE-BODY SCATTERING EQUATIONS 

We solve the Faddeev equations for three identical particles in the form 

T\(j>) =tP\(j>)+tG Q PT\(j)) (2.1) 

where t is the two-body t-matrix defined in the subsystem and the operator P is the sum 
of a cyclic and a anticyclic permutation of three objects. The initial channel state \<f>) is 
composed of a deuteron \(p d ) and the momentum eigenstate |q ) of the projectile nucleon. 
The free three-nucleon propagator is given by Gq = (E — H + ie) -1 , with E being the total 
center of mass (cm.) energy 
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E=^- q % + E d . (2.2) 

Here Ed is the binding energy of the two-body subsystem. Since we work below the three- 
particle breakup threshold we only need to consider the operator for elastic scattering 

U = PG 1 + PT. (2.3) 

From U one obtains the differential cross section for elastic scattering || as 
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(qVd|*/|qo¥>d) 
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(2.4) 



where |q'| = |qo|. The total cross section is either obtained by integrating over the angle 
variable 



°tot 



-i 

or via the optical theorem 



(go^Vdl^lqoVd) 
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(2.5) 



4?Ti 

a tot = -(2vr) 3 — Im( (q ^ d \U\q ^d) ) • (2.6) 



3<7, 
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In order to solve (|2.1| ) we introduce the standard Jacobi momenta p, the relative mo- 
mentum in the subsystem, and q the relative momentum of the spectator with respect to 
the subsystem. With \(f>) = |qoy?d) Eq. (|2.1|) reads 

(pq|T|qo^ d ) = (pq|tP|q ^d) + (pq|tG FT|q oV 9 d ). (2.7) 

The driving term of Eq ( |2.1| ) is given by 

(pq|tP|</>) = J d 3 q' d 3 p' d 3 q" d 3 p" (pq|t|p'q'> (pVI^IpV) (pV» 

= J d 3 q'd 3 p'd 3 P " ( P q|Vq'> <p'q'l^|p"qo> (p"M- (2.8) 

The momentum states are normalized according to (p'q'|pq) = S 3 (p' — p) £ 3 (q' — q). To 
evaluate the permutation operator P = P12-P23 + -P13-P23 explicitly, the Jacobi coordinates 
in the different subsystems (12) and (13) need to be expressed through those defined in the 
subsystems (23), which gives 

qi = -p 2 - |q2 

Pi = -5P2 + f qa 
qi = Ps - |qa 

Pi = -|P3 - |Q3. (2.9) 
Then the permutation operator occuring in Eq. fl2.8|) can be evaluated as 
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(p'q'liVV) = (P'q'|p"q">2 + (p'q'|p"q">3 

= $ 3 (p' + |q' + q") W - q' - iq") 

+ <* 3 (p' " |q' " q") W + q' + §q"), (2.10) 

The indices 2 and 3 indicate the corresponding subsystem (for more details see ||). Inserting 
this relation into Eq. reduces the driving term to 

(pq|t P\4>) = t s (p, |q + q ; E - £q 2 ) y> d (|q + §q |) ; (2.11) 

were t s (p, q, E) is the symmetrized two-nucleon t-matrix, 

t s (p, q, E) = t(p, q, £) + t(-p, q, £). (2.12) 

Since we neglect spin, the deuteron consists only of an S state, and thus the wave function 
depends only on the magnitude of the momenta. 

Carrying out a similar calculation for the integral term in Eq. ( |2.7| ) leads to the explicit 
form of the Faddeev equation 

(pq|T|q ^ d ) = ts(p, |q + q ; E - ^q 2 ) (f d (\q + |q |) 

+ / <rq - \ — - , q + ^q ,q T Wd . 2.13 

J E - ±(q 2 + qq" + q" 2 ) 

The transition operator T is needed for all values of q. Thus one encounters the pole of 
the two-body i-matrix at the bound-state energy E d . Extracting explicitly the residue by 
defining 

t B (p, q, E) = e _ Ea (2-14) 
and similarily for T, Eq. ( p.l3j ) can be written as 

(pq|T|q ^ d ) = 4(p, |q + qo;£- 4^g 2 )<^d(|q+ |qol) 
, ^Iq + q";^-^ 2 ) (qjJqVn^q^ (2 1 ^ 

+ i q E-±(? + q.<f + q>*) E-£q»*-E d + ie- l ' J 

This expression is the starting point for our numerical calculation of the transition amplitude. 
Correspondingly the operator for elastic scattering in Eq. ( |2.3|) reads 



(qVd|tf|qo¥>d) = 2 ^ d (±q' + q ) [E - ±(q' 2 + q'q + q 2 )) tp d (4 + l<lo) 

+ 2 J d 3 g'V d (iq' + q")(q , + iq",q"|T|q ^ d ). (2.16) 

The transition amplitude T as given in Eq. Q2.15| ) depends on the vector variables qo, 
q, and p. Going to cm. coordinates and choosing the z-axis in the direction of qo we are 
left with five independent variables. Those are the magnitudes of the vectors q and p, their 
angles with respect to the z-axis and the azimuthal angle (p pq between them. The vectors qo 
and q define the x-z plane in a cartesian coordinate system. In these variables the matrix 
element for the transition amplitude can be written as (p, x p , cos tp pq , q, x q \T\q Lp d ) with 
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V 

q 



IpI 
|q| 



x p = p-q 



cos 



cosv9(p,q) 



Pxy ' q 



xy ■ 



(2.17) 



The index xy denotes the projection of the vectors into the x-y plane. In order to obtain 
the matrix elements (pq|T|qoV3d) ; Eq. (|2.15 ) needs to be solved. For the integration we 
choose the z axis parallel to q. This implies that the azimuthal angle (p between (q + |q") 
and q" for (q + |q", q"|T|q </? d ) in the kernel of Eq. ( |2.15|) is zero, and thus cos ip pq = 1 
in Eq. ( |2.17| ). This also means that we only need to solve Eq. fl2.15| ) for cosLp pq = 1, or 
in other words, that p lies in the same plane that is spanned by q and q". From these 
considerations follows that we only have to solve Eq. (|2.15Q for (p, x p , 1, q, x q \T\q (p^), i.e. 
for four independent variables instead of five, as it could be assumed from the considerations 
proceeding Eq. ( 2.17|) . Thus, for our calculations we arrive from Eq. ( 2.17 ) with cosip pq = 1 
at the following two additional angle variables 



x 
x 



pq 



p • q XpXg 
q ■ q". 



1 Xp \ l 1 Xg COS ^Ppq XpXg ~\~ \ l 1 Xp \l 1 Xg 



(2.18) 



The angle variables x p and x q are already given in Eq. ( |2.17 ). The momenta occurring in 
Eq. ( |2.15|) are now given explicitly as 



l|q + qol 
|q+ |qo| 



\q 2 + qq x q + ql ||q + q"| 



\q 2 



qq"x" 



-J/2 



q l + qqoX q + \ql 



|q+|q"l 



f + qq"x"+ \q" 2 . 



(2.19) 



Expressions for the remaining angle variables can be found in Appendix Inserting all 
variables into Eq. (|2.15|) , the final expression for the transition amplitude reads 



(p,x p , l,q,x q \f\q Lp d ) 



(2.20) 



t s (p,Jiq+qqoX q + q 2 , 



c^QXpq ~f" QftXp 

\q 2 + qq x q + q 2 



e - i^q 2 ) M^q 2 + qqox q + \ql) 



oo 1 2tt 1 ii 

+ / dq"q" 2 [ dx" [ d V " t s (p, ^\q 2 + qq"x" + q" 2 , 2^ + ^ ; E - 



-q z + qq"x' + q" 



3 „2 



■hn 



q z ) 



q 2 + qq"x" + \q" 2 , 



qx q + ^q"y qo 



X 



qq"x 



, ~L,q",y qo \T\qoPd) 



E-Mq 2 + qq"x" + q 



n2\ 



E ~ ^q" 2 



E d + ie 



This is a three-dimensional integral equation in four variables, namely p, x p , q, and x q . The 
advantage of our choice of the coordinate system is that the free propagator has a relatively 
simple form, it depends only on the magnitude of momenta and one angle. Though in 
the present work we stay with our calculations below the breakup threshold, this particular 
form of the propagator will be the most suited form for considering the solution of Eq. (|2.20|) 
above. 
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The matrix elements of T provide input to the calculation of the matrix elements 
(qVdl^lqo^d) according to Eq. ( j2.16| ). In this integration we choose the z-axis parallel 
to q', so that there is no azimuthal angle between (q' + |q") and q". This specific choice 
ensures that we only need T as function of four variables. For the explicit representation of 
U the follwoing angle variables are are needed, together with the magnitude of q' and q" 



cosv 



x = q • q 

x" = q'.q" = costf" (2.21) 



and the momenta 
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l|q' + qo| = \j\q a + v'lox' + ql = qo\l^ + 



|q' + |qol = v q' 2 + q'qox' + \ql = W ^ + 



l|q' + q"| = v%' 2 + i ,( i" x " + 1" 2 



|q' + §q"l = \jq 12 + q'q" x " + \q" 2 - ( 2 - 22 ) 

The explicit expressions for the remaining angles are calculated in Appendix [B]. The final 
expression for the elastic scattering amplitude is then given by 

(qVd|^|qo<Ai) = (?osVd|tf1qoVd) 

m 



OO 1 2-7T 

+ 2 / dq"q" 2 E _ JL E — / ^'Vd( v / |% 2 + ?og'V' + ^ ^ 



4m- 



1 



x< v % + +" ^" 2 , ^=£fl±i£L 1- '/"■ .v|V'|q„ r M). (2.23) 

/g 2 + g g'V + ±g" 2 



From (go^Vdl^lqo^d) we obtain the differential cross section according to Eq. (j2.4|), and 
the total cross section via Eqs. ( |2.5|) and (fJT 



III. CALCULATION OF SCATTERING OBSERVABLES 

For our model calculations Yukawa interactions of Malfliet-Tjon [pi type are used, 



" (p ' p) = h- (w^WTA ~ (p'-pVJ • (31) 

We study two different types of pairwise forces, a purely attractive Yukawa interaction and 
a superposition of a short-ranged repulsive and a long-ranged attractive Yukawa interaction. 
It should be pointed out that we calculate the potentials as functions of vector momenta and 
thus define the interaction as a truly local force acting in all partial wave. The parameters 
are given in Table |, which also lists the corresponding deuteron binding energies. The 
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parameters are chosen such that the deuteron binding energy is close to the experimental 
one. With these interactions we first solve the Lippmann-Schwinger equation for the fully- 
off-shell two-nucleon t-matrix directly as function of the vector variables as described in 
detail in Ref. ||. The resulting t-matrix is then symmetrized to get t s (p f ,p,x; E — -£^q 2 )- 
We would like to point out that after having solved the Lippmann-Schwinger equation on 
Gaussian grids for p, p', and x, we solve the integral equation again to obtain the t- matrix 
at points x = ±1. Thus, when solving Eq. fl2.15| ), we do not have to extrapolate numerically 



to angle points x of t s (p',p, x\E — -^q 2 ), which can very well be located outside the upper 
or lower boundary of the Gaussian angle grid of the t-matrix. 

The fully off-shell t-matrix, t(p',p, x, E), is obtained for each fixed energy on a symmetric 
momentum grid, which is divided as (0,pi) U (pi,p max ). The intervals contain NP1 and NP2 
Gauss points, with typical values of NP1 = 40 and NP2 = 16 points. Typical values for the 
interval boundaries are p\ = 20 fm -1 and p max = 60 fin -1 . For the angular integration x 
32 Gauss points are sufficient. Since the momentum region which contributes to a solution 
of the two-body t-matrix is quite different from the region of importance in a three-body 
calculation, we map our solution for t s onto a momentum grid relevant for the three-body 
transition amplitude. This is done by applying the Lippmann-Schwinger equation repeat- 
edly. The t-matrix t s (p',p, x, e) is obtained at energies e = E — ^-q 2 , exactly at the q values 



needed in the three-body transition amplitude of Eq. ( |2.2C ). For extracting the residue of 
the two-body t- matrix, Eq. (|2.14|) , we represent t B as 



(E-E d )t.-> V\<p d )(<p d \V. (3.2) 



The Malfliet-Tjon type potentials support only an s-wave bound state, and thus Eq. ( |3.2| ) 
reads explicitly 

(E-E d ) £ s (p,p') 

1 ljdp"p" 2 V (p,p")<p d (p")\ \]dp"p" 2 V (p',p")Mp")\ , (3-3) 



were Vo is the I = component of the potential. 

In order to solve Eq. ( |2.20| ) we follow the iterative procedure outlined in Ref. [IT]. The 



method consists of first generating the Neumann series of Eq. (|2.20p and then summing up 



the series using the Pade method [E^rjj]. We typically need to sum 15-18 terms to obtain a 
converged result. This is not surprising, since due to the presence of the three-body bound 
state the Neumann series itself diverges. 

The g'-integration in Eq. ( |2.20| ) is cut off at a value of g max = 20 fm -1 . The integration 
interval is divided into two parts, (0, gi) U (qi,q maK ), in which we use Gaussian quadrature 
with NQ1 and NQ2 points, respectively. The value for q\ is chosen to be 5 fm -1 . Typical 
values for NQ1 and NQ2 are 24 and 16. For the distribution of quadrature points we use the 



maps given in Ref. ||15|| . The x" integration requires typically at least 18 integration points, 
while for the ip" integration 16 points are already sufficient. The p variable is also defined 
in an interval (0,pi) U (pi,£> m ax), where pi is chosen to be 7 fm" 1 and p max = 30 fm -1 . The 
two intervals contain NP1 and NP2 points, and we usually choose NP1 = NQ1 and NP2 = 
NQ2, respectively. 
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When solving Eq. ([2.20 ) we have to carry out two-dimensional and three-dimensional 



interpolations on t s and T. We use the Cubic Hermitean splines of Ref . [Tj| . The functional 
form of those splines is described in detail in Appendix B of this reference and shall not 
be repeated here. We find these splines very accurate in capturing the peak structure of 
the two-body t-matrix, which occurs for off-shell momenta p ~ p' . An additional advantage 
of the Cubic Hermitian splines is their computational speed, which is an important factor, 
since the integral in Eq. ( |2 . 2 U| ) requires a very large number of interpolations. 

In Fig. [I] the real and the imaginary parts of the scattering amplitude T(p, x p = 
l,cosip pq = l,q,x q = l,<7o) as obtained from MT-IVa potential are displayed. The pro- 
jectile energy is 3 MeV, and the amplitude is taken in forward direction, i.e. the two angles 
are set to zero. The figures show that most of the structure of the amplitude is concentrated 
at small momenta p and q. The corresponding amplitudes derived from the MT-IIIa poten- 
tial are shown in Fig. [2|. Though the imaginary part has a little more structure for p smaller 
than 1 fm -1 , the function is in general very smooth. The real real part of T have for both 
potentials a quite similar structure. 

The solution for the transition amplitude serves as basic input to obtain the elastic 
scattering amplitude according to Eq. ( |2.23| ). In carrying out the integrals we use the same 
grids as in the integral equation for T. The differential cross section obtained from the MT- 
IVa potentials is shown in Fig. |3] as function of the projectile laboratory energy E and the 
scattering angle d for energies from 0.01 Mev to 3.2 MeV. As expected, for very low energies 
the differential cross section is isotropic, which indicates that in a partial wave description 
only s-waves contribute. At about 1 MeV projectile energy the differential cross section 
starts to develop its more characteristic shape, namely in forward and backward direction 
and a minimum around fl = 100°. In Fig. ^ the differential cross section obtained from 
the MT-IIIa potential is displayed as function of the projectile energy and the scattering 
angle. The obvious difference with respect to Fig. ^| is that the magnitude of cr(d) obtained 
form the MT-IIIa is about 5 times larger than the one obtained from MT-IVa at small 
energies. The difference is related to the different values of the 3-body scattering length for 
the different potentials. The one for MT-IIIa turns out to be 2.034 fm and is much larger 
than the one for MT-IVa, which is 0.887 fm. These numbers are related to the different 3- 
body binding energies, which are -19.8625 MeV for the MT-IIIa potential and -25.1632 MeV 
for the MT-IVa potential. According to Ref. |T7[ the scattering length can be calculated via 



2tt 

a = —m (q ip d \U\q ^ d ) 



(3.4) 

90=0 



Due to the different scale the onset of a deviation from the isotropic cr($) is not so easily 
visible in Fig. |], but it also occurs at about 1 MeV. In order to better compare the differential 
cross section obtained from both potentials, this observable is shown in Fig. [| at 3 MeV 
for both potentials as function of the scattering angle. Here one can see that cr(i?) is larger 
especially in forward direction for the purely attractive potential MT-IVa. This behavior is 
even clearer visible in Fig. |6], where the total cross section atot obtained from MT-IIIa and 
MT-IVa are shown as function of energy. 

For demonstrating and discussing the numerical stability and accuracy of our algorithm 
we choose the MT-IVa potential and a fixed energy and discuss the behavior of the observ- 
ables as function of the grid points. Table |T| contains the total cross section calculated at 
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3 MeV using Eq. (|2.5|) and also via the optical theorem, Eq. ( |2.6| ). The first calculations 
listed in the table are performed with a moderate amount of grid points in all variables. We 
then successively in increase the points one variable at a time and see that the variation in 
the total cross section stays about 0.3%. From this we conclude that our algorithm is very 
stable, and the numerical error, which necessarily has to occur due to the large number of 
interpolations is certainly not higher than 1%. 

We choose to consider the total cross section for this stability study, since state-of-the-art 
measurements of the total cross section have an accuracy of about 0.5% [Q. We also see 
the different ways of calculating <7 tot , differ consistently by about 0.2%, almost independent 
of the number of grid points used. From this we can conclude, that the general error of 
our calculation is 0.5% or lower, which is for all practical purposes sufficiently accurate. A 
further test of the accuracy and convergence of our numerical calculation is in the insertion 
of our converged solution for the transition amplitude T a further time into Eq. (|2.20|) and 
then recalculate the observables with this new solution for T. The comparison is carried out 
for the differential cross section at 3 MeV and listed in Table |TTT| . Here we used the results 
of the calculation with the highest number of points from Table [D]> but results for the other 
calculations are similar. As can be seen the agreement of the two calculations is excellent, 
and we conclude that our calculations are properly converged. 



IV. A SECOND CHOICE OF VARIABLES FOR THE SCATTERING 

EQUATIONS 

In Section |I| we described our choice of the coordinate system used to solve the integral 
equation, Eq. ( |2.13| ), for the transition amplitude T. There we choose the z-axis for the 
integration parallel to q, which has the advantage of giving the free propagator in the kernel 
a relatively simple functional form. It also led to Eq. ( p. 20 ) being a three-dimensional integral 



equation in 4 variables. Obviously, the above described choice is not the only one. In order 
to test our calculations we additionally solve Eq. (|2.13 ) with a different choice of coordinate 



system. For this specific calculation we choose the z-axis for the overall coordinate system to 
be parallel to qo as well as for the integration in Eq. (|2.13| ). Due to the rotational invariance 
of the problem we can choose the azimuthal angle between those 2 coordinate systems to be 
zero. With these assumption the variables necessary to explicitly express Eq. ( |2.20| ) are 

P = IpI 
q = |q| 

x P = p- qo 
Xg = q • qo 

cos ip pq = cos <f(p, q) = p xy ■ q xy 



Uq — qq" = x q x " + \A ~ xiy/l — x" 2 cos ip" 



V P = pq" 

4' = |q"| 
X = q ■ q 



x p x" + \fl — x 2 \/l - x" 2 cos((p pq - if") 



cost/ = cosy?"(q",q). (4.1) 
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The calculation of the remaining angles and momenta is straightforward and similar to the 
ones given in Section Thus they are not given here. Using the above definitions of the 
variables we finally arrive at the explicit expression for the transition amplitude T 



(p,x p ,cosip pq ,q,x q \T\q ip d ) = (p d (^ q 2 + qq x q + \ql) 
X i(p, J\q 2 + qq x q + ql ^ + S ■ E - £q 2 ) 



\q 2 + qq Q x q + ql 



OO 1 2-7T | 

+ / dq"q" 2 / dx" J d<p» E _j_ (q2 + qql/ 



y q + q" 2 ) 



\qx vq + q"y p 



xWsJ^q +qq y q + q , , E - ^q ) 

'\q l + qq"y q + q" 



1 „ll 



E - £- t q" 2 -E d + ie 



(4.2) 



As in Eq. ( |2.20|) we solve for T, where the residue at the deuteron pole is explicitly taken 



into account as described in Eq. ( |2.14j ). In the form of Eq. ( |4.2| ) the three-body propagator 
has an explicit angle dependence. In addition one has a three-dimensional integral depending 
on five variables. The latter makes the numerical solution an order of magnitude more time 
consuming. Thus Eq. ( f4.2|) is solved on similar grids as Eq. ( [2.20D , however with fewer grid 
points. 

After solving for T, we obtain the elastic scattering amplitude by employing Eq. ( |2.16| ). 
Using the same coordinate system, namely the z axis being parallel to qo, the explicit 
expression for U reads 

(qVdl^iqo^d) = (qox'(p d \U\q (p d ) 
= 2rt(q JfT^)(E-£(2 + x>)) 

oo 1 2tt 

+ 2 J<W<r E _j_J_ E +%e J dx" J WMjiti + + 
o 4 ™ y d -1 o 

qox' + h"x" 



x ( V% 2 + + l = 2 = „ , cos <p } q\ x"\T\ WA ). (4.3) 

'ql + qog"y + k" 2 



Here x' — q' • qo and y — q'- q". The calculation of the azimuthal angle </> between (q+ |q") 
and q" is more complicated and given in Appendix 



It should be clear from the beginning, that the solutions of Eq. (4.2) is not only much more 
time consuming, but will also a priori contain a larger numerical error due to the increased 
number of interpolations. In addition four dimensional interpolations are required, whereas 
for the solution of Eq. ( |2.20| ) the maximum dimension for the interpolation is three. In Fig. 
we compare the differential cross sections at 3.0 Mev obtained from both algorithms using 



a medium number of grid points only. As an aside, solving Eq. (i.2) with the same high 



number of grid points as in Eq. (|2.20|) is too expensive, especially since we only had in mind 
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to perform a rough comparison of the two schemes. Thus, we also did not perform the same 
amount of accuracy tasks as described in Section 111 for the solution of Eq. ( j2.2C ). As seen 
in Fig. [7|, both solutions are reasonably close, and the accuracy is good enough to establish, 
that in general both methods give similar results. However, for practical calculations, the 
procedure described in this section should not be recommended. 



V. SUMMARY 

An alternative approach to the state-of-the-art three nucleon scattering calculations, 
which are based on solving the Faddeev equations in a partial wave basis, is to work directly 
with momentum vector variables. We formulate the three-body scattering equations below 
the three-body breakup threshold for identical particles as function of vector Jacobi momenta 
and the projectile momentum specifically the magnitudes of the momenta and the angles 
between them. We would like to point out, that our specific formulation and the choices 
of coordinate systems is also applicable above the breakup threshold. However, here the 
logarithmic singularities, inherent to the breakup, have to be treated explicitly. 

As two-body force we concentrate on a superposition of an attractive and a repulsive 
Yukawa interaction, which is typical for nuclear physics, as well as on an attractive Yukawa 
interaction. The corresponding two-body t-matrices, which enter the Faddeev equations was 
also calculated as function of vector momenta. We neglected spin degrees of freedom in all 
our calculations. 

In order to obtain scattering observables, which are in our case the differential and 
the total cross section, one solves first an integral equation for the transition amplitude 
T. The scattering amplitude is then obtained by an additional integration over the half- 
shell amplitude T. This set of equations contains in essence four vector momenta, the 
projectile momentum, the Jacobi momenta, and a momentum vector as integration variable 
in the kernel. In principle, one has different choices of the coordinate system, in which the 
calculations are carried out. We present different choices, one leading to a three-dimensional 
integral equation in four variables for the transition amplitude T and one leading to a three- 
dimensional integral equation in five variables for T. Obviously, the first choice is the 
preferred one. It has the additional advantage that the free three-body propagator acquires 
a relatively simple form, which will become relevant considering scattering above the breakup 
threshold. 

Using the transition amplitude given as function of four variables we calculate the observ- 
ables for different projectile energies and test the accuracy and stability of our algorithms. 
We establish that our calculations have an overall accuracy of less than 0.5%, which is suf- 
ficient for all practical purposes, i.e. comparison with experimental measurements. We also 
calculate the scattering observables at one energy using the transition amplitude given as 
function of five variables. The two different algorithms are in qualitative agreement, which 
gives us confidence, that our calculation is correct. 

Summarizing we can state that the Faddeev equations for scattering below the breakup 
threshold can be handled in a straight forward and numerically reliable fashion when using 
vector momenta as variables. Our formulation allows to treat the logarithmic singularities 
above the breakup threshold in a straight forward fashion, and work along this line is in 
progress. 
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APPENDIX A: VARIABLES FOR THE TRANSITION AMPLITUDE T 



In this case we need to rotate around the y = y" axis by the angle d q . The vectors needed 
in the new coordinate system are 



Q< 



!o = i2(t?,)qo 



= go 



/ — sin u q 


\ cos i? g 



Q" = R(# q )q" = q" 



J 



I sin •&" cos ip" \ 

sin sin ip" 
\ cos d" ) 



(Al) 



/ cos d q sin § p — sin d q cos d p \ 
P = R($ q )p = p 

\ sin d q sin $ p + cos "§ q cos $ p J 

It is then straightforward to obtain the following angles: 



(A2) 



x" 2 cos p" 



y p = P • Q" = x"x pq + Vl - x" 2 (x q ^Jl - x 2 p - x p ^Jl - x 2 ) cosp" 
y qo = Q" ' Qo = x q x" - j\-x* q y/l = 

P . (IQ?Q") = h^ + ^'yp 

^\q 2 + qq"x" + q" 2 
Qo • (Q +]Q") - qXq + ~ 2q " Vqo 



q 2 + qq"x" + \q" 2 



(A3) 



APPENDIX B: ANGLES FOR THE SCATTERING AMPLITUDE U 

We need the angle between q and q". It can be calculated in terms of the integration 
variables in the new coordinate system. We have chosen the z axis for the integration being 
in the x-z plane of the original coordinate system, and the azimuthal angle between them 
is zero. To get to the new coordinate system we only need to rotate around y = y" axis by 
the angle d 1 using the rotation matrix 



/cos??' -sintf'N 

1 
\ sin §' cos J 



(Bl) 
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The two vectors in the new coordinate system are given by 



Qo = R(f)qa = q 



/ — sin i?' 


y cos •& 



Q" = RW = q" 



I sin d" cos ip" 
sin d" sin p" 
\ cos d" 



The angle between them is 

y = Q" ■ Qo = - sin sin d" cos </ + cos cos # 



x'x" — V 1 — X /2 Vi — X"^ COS C/9" 



o"2 ™ 



The other angle we are looking for is 



Qo-(Q' + |Q") 



q x' + \q"y 



ql + g g'V + \q" 2 



fB2l 



(B3) 



(B4) 



APPENDIX C: ANGLES IN THE FIVE DIMENSIONAL CASE 



We are looking for the azimuthal angle p> between (q + |q") and q". We had chosen q 
being in the x-z plane, now one can find the azimuthal angle p> by going into the x-y plane 
and using the cosine theorem 



h"J 2 + l(q + k'W 2 



2||<4I l(q+!q")J 



I Qxy | i 2 "-ixy 

The components of the vectors in the x-y plane are given by 



cos p. 



I *-lxy 1 

-a" I 

1„/A 



gsin-iL = q Jl — x 



(q + V) 



ig"sin^ g » = — x 

|q+|q"|sin^ q0i( i q „ +q) 



//2 



|q+k"| 



1-x? 



with 



a?i = qo ■ (q + jq") 



\q"x" + qx g 
|q+lq"| ■ 



Inserting these relations in Eq. ( |CT| ) we find 

^ L |q+|q'f (i-x\)-q 2 (i 



cos p 



\q"\l 



x 



x„ 



q"Vl ~ x //2 |q+ iq /, |v / T^i' 
which can be simplified to the final result 



1 q'Wl - x" 2 + qJl - x 2 cos p" 



cos p 



X' 



n2\2 



2^2 



+ (g'Vl - x" 2 )(qJl - x 2 ) cos ip" 



(CI) 



(C2) 



(C3) 



(C4) 



(C5) 
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TABLES 



TABLE I. Parameters and deuteron binding energy for the Malfliet-Tjon type potentials. 
As conversion factor we use units such that He = 197.3286 MeV fm = 1. We also use 



h 2 /m = 41.47 MeV fm 2 . 





V A [MeV fm] 


HA [fur 1 ] 


Vr [MeV fm] 


HR [fm x ] 


E d [MeV] 


MT-IIIa 


-626.8932 


1.550 


1438.7228 


3.11 


-2.231 


MT-IVa 


-65.1776 


0.633 






-2.223 



TABLE II. The total cross section for the MT-IVa potential at 3.0 MeV. The number of grid 
points for f are NQ = NQ1 + NQ2, NP = NP1 + NP2, NXP, and NXQ as explained in the text. 
The grid points for the integration are NX" and Nip". The grids for the t-matrix are denoted with 
NP t and NX t . 



NQ 


NP 


NXP 


NXQ 


NX" 


Nip" 


NP t 


NX t 


int. 
°tot 


opt. 
^tot 


30 


30 


18 


14 


18 


10 


40 


40 


3903.18 


3906.53 


30 


30 


18 


14 


18 


16 


40 


40 


3909.12 


3913.19 


30 


30 


18 


14 


18 


20 


40 


40 


3909.22 


3913.95 


30 


30 


18 


18 


18 


20 


40 


40 


3090.93 


3914.09 


30 


30 


18 


22 


18 


20 


40 


40 


3910.06 


3914.09 


30 


30 


22 


22 


18 


20 


40 


40 


3909.73 


3914.00 


30 


38 


22 


22 


18 


20 


40 


40 


3909.66 


3914.30 


30 


42 


22 


22 


18 


20 


40 


40 


3909.11 


3913.45 


38 


42 


22 


22 


18 


20 


40 


40 


3916.66 


3921.11 


42 


42 


22 


22 


18 


20 


40 


40 


3913.17 


3917.47 


46 


42 


22 


22 


18 


20 


48 


40 


3911.16 


3915.16 


46 


42 


22 


22 


18 


20 


56 


40 


3911.13 


3915.02 


46 


42 


22 


22 


18 


20 


64 


40 


3911.11 


3914.96 


50 


42 


22 


22 


18 


20 


64 


40 


3911.15 


3915.06 


50 


42 


22 


22 


24 


20 


64 


40 


3910.88 


3915.04 


50 


50 


22 


22 


24 


20 


64 


40 


3910.78 


3914.93 


50 


50 


22 


22 


24 


22 


72 


40 


3910.79 


3914.88 


50 


50 


22 


22 


24 


22 


80 


40 


3910.71 


3914.72 


50 


50 


22 


22 


24 


22 


80 


40 


3910.71 


3914.69 


58 


58 


22 


22 


24 


22 


80 


40 


3913.31 


3917.22 


58 


58 


26 


26 


24 


26 


80 


40 


3913.95 


3917.44 


58 


58 


30 


30 


24 


26 


80 


40 


3913.85 


3917,44 


58 


58 


34 


34 


24 


26 


80 


40 


3913.85 


3917.46 


58 


58 


38 


38 


24 


26 


80 


40 


3913.86 


3917.61 
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TABLE III. The differential cross section for the MT-IVa potential at 3.0 McV. The second 



column contains the cross section obtained from the Pade sum and the third column the cross 
section obtained by reinserting the solution. The corresponding total cross sections are 3913.86 
mb and 3913.93 mb respectively. 



[deg] 


o-(tf) [mb] 


ff(tf) [mb] 


0.000 


634.458 


634.483 


8.215 


624.423 


624.446 


17.548 


590.018 


590.040 


26.893 


535.377 


535.396 


36.240 


467.039 


467.056 


45.589 


392.139 


392.153 


54.938 


316.841 


316.853 


64.288 


245.326 


245.336 


73.638 


180.241 


180.249 


82.988 


123.550 


123.557 


92.337 


78.909 


78.917 


101.687 


53.890 


53.899 


111 H'J'7 

111.03 7 


61.6 to 


61.690 


120.387 


120.744 


120.760 


129.736 


250.915 


250.937 


139.085 


464.483 


464.511 


148.434 


752.825 


752.860 


157.780 


1074.342 


1074.383 


167.122 


1354.729 


1354.775 


176.421 


1508.989 


1509.036 


180.000 


1522.727 


1522.774 
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FIGURES 



Re(T) [10~ 9 fm 3 ] 




FIG. 1. Real and imaginary part of T(p, x p = 1, cos (p pg = 1, q, x q = 1, qo) at 3 MeV projectile 
energy as obtained from the MT-IVa potential. 
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Re(T) [10~ 9 fm 3 ] 





FIG. 2. Same as Fig. 1 but for the MT-IIIa potential. 
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FIG. 3. The differential cross section <j($) as function of the projectile energy and scattering 
angle $ obtained form the MT-IVa potential. 




FIG. 4. Same as Fig. 3 but for the MT-IIIa potential. 
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FIG. 5. The differential cross section a(-&) at 3 MeV projectile energy obtained from the 
MT-IIIa potential (solid line) and the MT-IVa potential (dashed line). 
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FIG. 6. The total cross section <Jtot as function of the projectile energy obtained from the 
MT-IIIa potential (solid line) and the MT-IVa potential (dashed line). 
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20 40 60 80 100 120 140 160 180 

FIG. 7. The differential cross section obtained from the MT-IVa potential at 3.0 MeV based 
on the solution of T using four variables (solid line) and using five variables (dashed line). The 
explanation of the calculations is given in the text. 
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